Spin Transport in the XXZ Chain at Finite Temperature and Momentum 
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We investigate the role of momentum for the transport of magnetization in the spin-1/2 Heisen- 
berg chain above the isotropic point at finite temperature and momentum. Using numerical and 
analytical approaches, we analyze the autocorrelations of density and current and observe a finite 
region of the Brillouin zone with diffusive dynamics below a cut-off momentum, and a diffusion 
constant independent of momentum and time, which scales inversely with anisotropy. Lowering the 
temperature over a wide range, starting from infinity, the diffusion constant is found to increase 
strongly while the cut-off momentum for diffusion decreases. Above the cut-off momentum diffusion 
breaks down completely. 
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Understanding spin transport in quantum many- 
particle systems is a fundamental challenge to physics, of 
immediate relevance to future information technologies 
[l|, and intimately related to timely issues of dynamics 
and thermalization in a more broader context [2j. While 
conventional spin conductors like silicon [3|, III-V semi- 
conductors [J] , carbon nanotubcs [5| , or graphene [6| nec- 
essarily feature spins which are associated with itinerant 
charge carriers, insulating quantum magnets may open 
new perspectives for spin transport, with pure magneti- 
zation currents flowing solely by virtue of exchange in- 
teractions. Magnetic transport in one-dimensional (ID) 
quantum magnets has experienced an upsurge of inter- 
est in the last decade due to the discovery of very large 
magnetic heat conduction [7| with mean free paths above 
l/im [8|. Genuine spin transport in quantum magnets 
remains yet to be observed experimentally, however long 
nuclear magnetic relaxation times [9( have been estab- 
lished, which even allow for manipulation with magnetic 
fields 0. 

Theoretically, significant attention has been devoted to 
spin transport in ID quantum magnets, see Refs. [ll| and 
Il2i for reviews. The dissipation of spin currents is a key 
issue in this context and has been analyzed extensively 
at zero momentum and frequency in connection with the 
spin Drude weight [13| . Spin current dynamics at finite 
momentum remains one of the open questions. In this 
Letter, we will address this question for the antiferromag- 
netic and anisotropic spin-1/2 Heiscnbcrg (XXZ) chain 

N 
r 

where S l r (i = x,y,z) are the components of spin-1/2 
operators at site r, N denotes the number of sites, J > 
represents the exchange coupling constant, and A is the 
anisotropy. The XXZ chain is a fundamental model to 
describe magnetic properties of interacting electrons. It 
is relevant to the physics of low-dimensional quantum 
magnets [141 ], ultra-cold atoms [15|, nanostructures f 161 ] - 



and - seemingly unrelated - fields such as string theory 
[17| and quantum Hall systems [18| ■ 

Early analysis of the time-dependent correlation func- 
tion of the local spin density has been performed in the 
high-temperature limit, T — oo, sug gesting the absence 
of spin diffusion for ^ A ^ 1 [19j. Subsequent, stud- 
ies have concentrated on the spin Drude weight at zero 
momentum q = [131 ] , allowing for no conclusions on dif- 
fusion laws at finite momentum. First low-temperature 
quantum Monte-Carlo studies at q ^ [20( found no ev- 
idence for spin diffusion; however, more recent results 
from bosonization and transfer-matrix renormalization 
group [21| as well as quantum Monte-Carlo |22j are con- 
sistent with finite-frequency spin diffusion in the small- 
momentum regime, at A = 1 and for low temperatures 
T <C J , with a spin diffusion constant D which diverges 
oc 1/TlnT. The physics at intermediate temperatures 
and arbitrary momenta remains undisclosed. 

Therefore, in this Letter, we consider the transport of 
magnetization by analyzing autocorrelations of spin den- 
sity and current at finite momenta, covering the com- 
plete Brillouin zone, and at intermediate temperatures 
0.5J ^ T < oo (h = ks = 1). We focus on the case of 
finite anisotropy A > 1, where Eq. (fT]) features a gapped 
ground state. Using a combination of exact diagonal- 
ization and perturbation theory, we uncover a regime of 
diffusive transport below a finite critical momentum qr>. 
In this regime, density modes at fixed momentum q decay 
with a diffusion constant D q and our analysis is consis- 
tent with D q independent of momentum and inversely 
proportional to the anisotropy. As the temperature is 
lowered from T = oo, we observe a decrease of the criti- 
cal momentum and an almost exponential increase of the 
diffusion constant. We provide evidence for a complete 
breakdown of diffusion above the critical momentum. 

We begin by introducing the generalized diffusion co- 
efficient as a quantity suitable to describe the evolu- 
tion of a harmonic spin density profile close to equilib- 
rium, i.e., in the linear response regime. To this end, 
the central quantities we analyze are the autocorrela- 



tion functions C s . q {t) = Rc(S*(t)Si q )/N and Cj, q (t) = 
Rc(J^{t)Jl q )/N of the spin density S z q = £ r e l i r S* and 
the spin current J* = JJ\ e % i r (S*S v r+l - SyS* +1 ) at 
momentum q = 2irk/N [231 ] . where Re indicates the 
real part, (...) denotes the canonical equilibrium aver- 
age at the inverse temperature j3 = 1/T, and t rep- 
resents the time. Since the density S q and the cur- 
rent J q are connected by the lattice continuity equa- 
(1 — e lq )J q , the autocorrelations are re- 
-q 2 C, 



1.2 



tion dtSq = 

lated by d 2 Cs, q (t) = —q 2 Cj^ q {t) with the abbreviation 
q 2 = 2(1 — cosg). The generalized, time- and momentum- 
dependent diffusion coefficient is defined via 

d t C s , q (t) W 



D q (t) 



-q 2 C, 



S,S) Cs, q (0) - q 2 I 2 q {t) 



(2) 



To 



arrive at 
we integrate 
dtC s , q (t)\ t 



=o 



the right-hand expression in Eq. ([2"]l. 
the continuity equation twice, using 
= and introducing the two integrals 
I\(t) = f Q &t>Cj, q {t>) and I 2 (t) = H&t'Ilit'). 

The left-hand expression in Eq. © identifies the quan- 
tity q 2 D q (t) with the instantaneous decay rate, at time 
t, of a spin density profile with wave vector q close 
to equilibrium. Fick's law corresponds to the case of 
D q (t) = const. The main goal of this Letter is to ana- 
lyze the time- and momentum-dependence of this quan- 
tity versus temperature. We emphasize that a complete 
knowledge of this dependence allows to propagate arbi- 
trarily shaped spin density profiles in time. This does 
not only share a common interest with time-dependent 
density-matrix renormalization group studies [24[, yet 
confined to zero temperature, but even more so may be 
of relevance to laser pulse induced time-dependent trans- 
port measurements, including recently proposed time-of- 
flight and thermal imaging techniques [251 ] . 

Qualitatively, the variation of D q (t) versus t can be 
understood from a standard relaxation-time approxima- 
tion, in which the current autocorrelation Cj yq (t) = 
cxp(— t/t q ) Cx g (0) decays exponentially. For short times, 
t < t q , Eq. © then yields D q (t) - 1 - e - */**, which 
starts with a linear increase, D q (t) oc t, and turns into a 
'plateau' D q {t) « const., starting at t — r q > t q . This 
plateau marks the hydrodynamic regime. Namely, pro- 
ceeding to the long-time limit, i.e. for t 3> t q , and to the 
long-wavelength limit, i.e. for q 2 (t - t q )D q < 1, Eq. @ 
leads to a time- independent diffusion constant D q (t) = 

0)/C s , g (0), which is 



tqCj, 



D q=0 . 



D a + 0{q 2 ), where D q 

equivalent to Einstein's relation [261 ] . and Do 

In principle, partial conservation of currents at q = 0, 

i.e. the impact of a finite Drudc weight at zero frequency 

[131 ] . can also be included into this qualitative picture. 

For that case the exponential decay of Cj${t) has to be 

leveled off into Cjfl(t — > oo) = const. > 0. This leads to 

a linear increase Do(t — > oo) oc t. However, the Drudc 

weight will not be an issue in this Letter. In fact, there 

is no zero- frequency contribution of currents at q ^ 0, 

which follows directly from the continuity equation. 
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FIG. 1. (color online) The time- and momentum-dependent 
diffusion coefficient D q (t) at j3 — and (a) A = 1.5, 
(b) A = 2.0. ED results are shown for N = 18 and 
q/(2ir/N) = 0,1,2,3,4, and 9 (non-solid curves). POPT re- 
sults are shown for q = (red/dark-colored, solid curves) and 
q > (green/light-colored, solid curves). Thick arrows on the 
ij-axis mark the locations of the current decay time r q . 



While the gross feature of the preceding relaxation- 
time ansatz can serve as a guideline to interpret the re- 
sults of unbiased exact diagonalization data, on which we 
will report later, it is not justified a priori. Therefore, 
and to gain a deeper insight into the high-temperature 
current dynamics generated by the Hciscnberg model, 
we will first turn to a quantitative discussion using an 
analytical method. This method employs the projec- 
tion operator perturbation theory (POPT) of Ref. |27I . 
which allows to derive a rate equation d t Cj[s], q {t) = 
—lJ[s].q{t)Cj[s], q {t) for the current [density] autocorre- 
lation. This rate equation gives access to D q {t) through 
the right-hand [central] expression in Eq. ©. The POPT 
yields a short-time expansion for the decay rate ~/j[s] tg {t), 
the terms of which can be evaluated from a decomposi- 
tion H = Ho + Hi, if the observable of the autocorre- 
lation Cj[s], q (t) is a conserved quantity for the unper- 
turbed Hamiltonian Ho. For the current, we choose the 
XY model for Ho, in which J q is conserved only at q = 0. 
For the density, we choose the Ising model for Ho, in 
which S q is conserved for all q. Then for short times we 
obtain approximately: 
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For the full quantitative evaluation of D q (t) we determine 
the leading-order term in Eqs. ( El an d (j4]) numerically ex- 
act, following the scheme in Ref. [27|, which leads to small 
changes only. We note that for a complete integration of 




FIG. 2. (color online) Finite-size scaling results for the dif- 
fusion coefficient D q (t) at (a) q — 0, (b) q/(2n/N) = 1 for 
different N = 10, 12, . . ., 18 at /3 = and A = 1.5. In (a) 
finite-size variations can be neglected for t J < 10 at N — 18. 
While q — const, cannot be maintained in (b), the tendency 
is similar to (a). 



Eq. ([2]) the high-temperature limits of the static correla- 
tion functions are needed, i.e. CVgi g (0) = 1/8 [1/4]. 

For q — 0, we obtain from the POPT and the right- 
hand expression in Eq. ([5]) a leading order prediction as 
follows: the current autocorrelation Cjfi{t) decays even 
stronger than in a simple relaxation-time approximation, 
i.e. according to a Gaussian, and the diffusion coefficient 
Do(t) is an error function. This prediction is consistent 
with using Eq. ([3]) for all times, which is justified because 
the current relaxation time from Eq. ([3]) is to J w 1.9/ A, 
C*j,o(io)/Cj,o(0) = 1/e. Therefore for A = 1.5 or 2.0, 
as in this Letter, Do(t) has saturated for times within 
tJ < 1.5. The resulting quantitative Do(t) is shown in 
Fig. Q] (red/dark-colored, solid curves): Here, D (t) first 
increases linearly but then saturates at a constant value 
D Q J « 0.88/A, which is reached at tJ > t q J rj 3.0/A. 
For the remainder of this Letter we refer to the satura- 
tion time T q as the 'current relaxation time', rather than 
t q , since it can be extracted more precisely from later nu- 
merical data. We emphasize that our value of D p ag rees 
remarkably well with other approaches in Refs. [28| and 
29i . It is worth to mention that, for A — > oo, the 1/A- 
scaling of Dq may break down due to possible recurrences 
of Cj.o (t), see also Ref . [32| for an alternative point of view. 

For q ^ 0, we obtain from the POPT and the central 
expression in Eq. (J2)) a prediction for the full momentum- 
dependence of D q (t). This prediction is only valid at 
short times, set by Eq. (QJ. The resulting quantita- 
tive D q (t) is depicted in Fig. Q] (green/light-colored, solid 
curves). Clearly, D q (t) is not constant in the short-time 
domain as a function of q. The g-dependence arises from 
the next-to-leading order term of the POPT and becomes 
significant for momenta above q ~ 0.27rA, at A = 1.5 
and 2.0, and is particularly evident for q = n. 

In order to complete the picture at j3 = for arbi- 
trary momenta and times we apply exact diagonalization 
(ED) to chains of length N = 18, allowing for a g-grid 
with Sq w 0.117T. Figure Q] depicts our results for D q (t). 
Several comments are in order. First, Figs. [TJ (a) and 
(b) show a convincing agreement between ED and both 



POPTs within their respective ranges of validity, which 
corroborates our analysis. Next, a given density mode 
at wave vector q shows the signature of a diffusive de- 
cay if there is a plateau with D q (t) rj const, within a 
'long-time' window r q < t < tjj with r q <C trj. Clearly, 
Fig. Q] shows that the first three (four) momenta for 
A = 1.5 (2.0) feature such plateaus. Long-time devia- 
tions from this behavior can have several origins, such 
as finite-size effects, finite Drudc weights, or other low- 
frequency anomalies. Most remarkable, the plateau val- 
ues of D q (r q < t < to) ~ Dq and r q m tq are inde- 
pendent of momentum, to within the typical finite-size 
variations which occur for N = 16 — > 18. Finite-size 
effects are further quantified in Fig. [5] for the interval 
tJ < 10. While a constant finite q cannot be main- 
tained as N varies, it is still obvious that system sizes 
of N = 18 are completely sufficient to determine D q 
at the plateau. The weak dependence on q has to be 
contrasted against the significant q-dependence for larger 
momenta. In agreement with the POPT for q = 0, the 
values of Do and To from ED can be scaled onto a sin- 
gle expression for the two anisotropics studied, namely, 
D Q J « 0.88/A for t J > t J « 3.0/A. This is one of 
the main results of this Letter, i.e., the existence of an 
extended momentum-space region with a q- independent 
diffusion constant oc 1/A. Clearly, the number of mo- 
menta to which the diffusion criterion applies is smaller 
in Fig. [1] (a) than in (b). For both A = 1.5 and 2.0, we 
find no indications of diffusion for q > 0.227rA = qrj. In- 
stead, D q>qD (t) displays divergent behavior due to oscil- 
lations of Cs, q {t) with time, preventing diffusive behavior 
to occur. These oscillations have already been reported 
in Ref. [30j for smaller A. 

We emphasize that ED results for the spectra Cj\s], q (w) 
at small q versus frequency uj agree with our interpreta- 
tion from the time domain. E.g., focusing on A = 1.5, 
Fig. [3] (a) shows that the spectrum Cj tq (u))/Cs, q (t — 0) 
is consistent with a Gaussian of height DqJ rj 0.59, as 
predicted by the POPT at q = 0. The low-frequency 
behavior is still governed by finite-size effects and de- 
viations from the Gaussian occur at a frequency scale 
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FIG. 3. (color online) Spectrum of the (a) current and (b) 
density autocorrelation at f3 — and A = 1.5. ED results are 
shown for N = 18 (symbols). In (a) a Gaussian with height 
Do and in (b) a Lorentzian with width q 2 Do are indicated for 
comparison (curves), using DoJ ~ 0.59 from the POPT. 
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FIG. 4. (color online) (a)-(c) ED results for the diffusion 
coefficient D q {t) at different /3 > for A = 1.5 and JV = 18 
(curves). Thick arrows mark the approximate locations of the 
current decay time r q . (d) The resulting diffusion constant Do 
versus j3 for A = 1.5 (squares) and A = 2.0 (circles), which 
can be determined with a precision of 20% (error bars). 



u>/J < 1, which is independent of q. This agrees with 
the g-indepcndent time scale in Fig. f2J where finite-size 
effects set in. Similar spectra of Cj^ q (uj) have been ob- 
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0. Note that the (finite size) 
= is not shown in Fig. [3] (a). 



tained in Ref. 

q = Drudc weight at ui 

Figure [3] (b) shows that Cs, q (u)) is consistent with a 

Lorentzian of width (J 2 Dq, again DqJ « 0.59, as expected 

for diffusive density decay. 

Now we turn to the effects of temperature by increasing 
j3 from to (3 J = 2. Since the POPT is not applicable at 
f3 y^ 0, we focus on the ED results. Figure E](a)-(c) sum- 
marizes our findings for A = 1.5 and (3 J = 0.5, 1, 2. We 
observe two effects. First, as the temperature is lowered, 
the number of momenta with diffusive density dynamics 
decreases. At (3J = 0.5 and 1 the mode with q — 0.117T 
still decays diffusively but for (3 J = 2 only the q = 
mode displays diffusion. Second, as the temperature is 
lowered, D q and r q increase significantly. For q = this 
increase can be followed up to (3 J = 2. Fig. @] (d) dis- 
plays Dq versus (3 in a semi- logarithmic plot for A = 1.5 
and 2.0. From this plot, one might be tempted to spec- 
ulate on an exponential increase of Dq with f3 beyond 
the temperature window depicted, see a related claim in 
Ref. l32i However, in view of the hydrodynamic relation 
Dq = t o Cjfi(0)/Cs,o(0) this is a subtle issue. From our 
numerical analysis, we find Cs,o(0) to be the dominant 
source of D 's T-dependence for < (3 J < 2. But Cs,o(0) 
is not tx exp(c/3) for all j3 [33J. An exponential increase of 
Do must further break down as j3 — > oo due to the finite 
spin gap for A > 1. We also mention that, for A = 1 and 
f3J 3> 1, the dominant T-dependence of Do stems from 
to oc l/(TlnT) [2lj,|22|, which is not exponential. 

Finally, we turn to a more detailed discussion of the 
temperature dependence of the critical momentum qo- 



(b) A=2.0 




FIG. 5. (color online) The diffusive range of momenta q versus 
f] for (a) A = 1.5 and (b) A = 2.0. ED results are shown 
for N — 16 (squares) and N = 18 (circles). [Open symbols 
are borderline values, see Fig. [JJ and Fig. U (a)-(c).] The 
curve q 2 Do to = 1 is indicated for comparison (solid curves) . 
(Dashed curves estimate errors for Do and to.) 



To this end, we first collect all momenta q < qn in Fig. [5] 
Then, to rationalize this, we invoke the standard hydro- 
dynamics criterion that the relaxation time l/(q 2 D q ) of 
a diffusive density mode should be larger than the decay 
time T q of the current, or equivalcntly, that a diffusive 
density spectrum should be narrower than the current 
spectrum, see Fig. |3j Therefore, breakdown of diffusion 
occurs at q 2 D q r q ~ 1, where we may set D q = Dq and 
T q = To , due to the weak g-dependence of these quantities 
in our case. Based on our ED results for D and t , Fig. [5] 
displays the lines q 2 D tq = 1 versus j3 for both A = 1.5 
and 2.0 (solid curves). The obvious agreement between 
these lines and the boundaries for the collected values of 
q < qo is a convincing consistency check of our approach. 
Apparently, as j3 increases, qD decreases. In view of the 
temperature dependence of Dq and To, this decrease is 
also approximately exponential for < f3J < 2. To asses 
the relevance of finite size effects, Fig. [5] contains a com- 
parison between the lines q 2 DoTo = 1 and the observed 
diffusive modes for N = 16 and 18 (symbols). Given the 
limited resolution of the g-grid, the agreement with these 
two system sizes is remarkably good. 

In summary we have investigated magnetization trans- 
port in the spin- 1/2 XXZ chain above the isotropic point 
at finite temperature and momentum. We found an ex- 
tended momentum-space region of spin-diffusion with an 
approximately time and momentum independent diffu- 
sion constant. The diffusion cut-off wave vector (diffu- 
sion constant) was found to scale approximately linear 
with the (inverse) anisotropy and to decrease (increase) 
strongly with the inverse temperature. 
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